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ABSTRACT 

We reanalyze the galaxy-mass correlation function measured by the Sloan Digital Sky Survey to obtain host 
dark matter halo masses at galaxy and galaxy group scales. We extend the data to galaxy clusters in the 
2MASS catalog and study the relation between central galaxy luminosity and halo mass. While the central 
galaxy luminosity scales as ^ ^0-7-o.8 masses, the relation flattens to ^ above ^ 4 x IO'^Mq. The 

total luminosity of galaxies in the halo, however, continues to grow as a power-law ^ ^> 8-0.9 starting from 
the hypothesis that the central galaxies grow by merging (“galactic cannibalism”), we develop a simple model 
for the evolution of their luminosities as a consequence of the accretion of satellite galaxies. The luminosity- 
mass relation flattens when the time scale on which dynamical friction induces orbital decay in the satellite 
galaxies exceeds the age of the dark matter halo. Then, the growth of the central galaxy is suppressed as it can 
cannibalize only the rare, massive satellite galaxies. The model takes the dependence of the total luminosity of 
galaxies in a halo on its mass and the global galaxy luminosity function as input, and reproduces the observed 
central galaxy luminosity-mass relation over three decades in halo mass, (10'^- 1O'^)M0. The success of 
the model suggests that gas cooling and subsequent star formation did not play an important role in the final 
assembly of central galaxies from sub-L* precursors. 

Subject headings: cosmology: observations — cosmology: theory — galaxies: clusters: general — galaxies: 
formation — galaxies: fundamental parameters 


1. INTRODUCTION 

Our understanding of galaxy formation and evolution is 
still incomplete. The mass function of dark matter halos 
is routinely measured in numerical simulations. If a linear 
relation is assumed between the galaxy luminosity and the 
halo mass, the abundance of galaxies at the faint and the 
bright en ds of the luminosity range is signiflcantly below the 
expec ted JVale & Ostrikeil2004lvan den Bosch, Yang, & Mo| 
20041). The standard picture of galaxy formation (e.g., 
Rees & Ostrikeill977HWhite & Reeal97ll) postulates that as 
dark matter halos virialize, the gas inside the halos is heated to 
the virial temperature and then cools to form galaxies. While 
the dearth of faint dwarf galaxies can be explained in semi- 
analytic models by invoking a variety of feedback mecha¬ 
nisms that expel gas from small dark matter halos at early 
times, these mechanisms eventually yield an overabundance 
of bright galaxies as the expelled gas event ually cools in 
massi ve halos to form luminous galaxies (e.g.. lBenson et all 

iinn). 

Motivated by the discovery that the temperature distribu¬ 
tion of g as in virializing halos in numerical simulations is 
bimodal JKatz et alJl2003^ . and by the lack of evidence for 
significant coo l ing flo ws in galaxy clusters ('e.g.. lEahian et alJ 
1200 Ih . iBinnevI j2004l) suggested that only the colder compo¬ 
nent cools to form galaxies, while the hotter component re¬ 
mains at the virial temperature. One mass scale characterizing 
galaxy formation would then be that below which the shocks 
cannot be sustained, the gas is not heated, and all the gas in 
the colder component streams to the halo center and forms a 
galaxy. Models of shock formatio n suggest this critical mas s 
is small, Mshock ~ (1“6) x 10 "M^ J Pekri & Birnboirr J2004l) . 
A similar result is derived in iMaller & BullockI ( |2004[) based 
on models of multi-phase galaxy formation. 

Since galaxies do not continue to grow by accreting hot 
gas, their stellar masses can increase only by accreting other 


galaxies. Thus, we expect that the galaxy growth on scales 
above Mshock is simply governed by the physics of dissipation¬ 
less merging of pre-existing galaxies. The merging of galax¬ 
ies tra cks the hierarchical assembly of the host dark matter 
halos llOs triker & Tremaine! 1975HOstriker & Hausmanll977t 

iMerrittUlQfl^TDubinsklTTQQflh . Here, a second characteris¬ 

tic scale is expected, which is associated with the efficiency 
by which galaxies accrete other galaxies. Above this second 
scale, smaller galaxies that have merged into a larger halo do 
not have enough time to reach the halo center and aggregate 
into one final central galaxy. 

In §|2j we extract the dependence of central galaxy luminos¬ 
ity on halo mass in galaxies imaged with SDSS and 2MASS 
surveys. We provide evidence for a characteristic scale in this 
dependence at ~ (1-6) x IO^^Mq. In §□ we derive a relation 
between the luminosity of the central galaxy and the mass of 
its host dark matter halo from first principles. The derived 
relation fits the data from galaxy scales ~ 10 '^Mq to cluster 
scales ^ IO'^Mq. 

We adopt the curre nt concordance cosm ological model con¬ 
sistent with WMAP dSnergel et al.ll2flfl3h . Unless otherwise 
noted, M refers to the virial mass of a halo out to an overden¬ 
sity of 200, and rvL is the associated virial radius. 

2. THE GALAXY-DARK MATTER CORRELATION EUNCTION 

The luminosities of central galaxies in clusters and groups 
can be extracted directly from IR imaging data in the 2MASS 
survey, wh ile the masses can b e inferred form the X-ray tem¬ 
peratures JLin &Mohdl2()()l . As a direct mass measure¬ 
ment of the halo mass is impossible at galaxy scales, we ex¬ 
tract the masses from tangential shear measurements in weak 
gravitatio nal lensing around a sample of foreground galax¬ 
ies (e .g., iBrainerd, Blandford. & SmaOl ll 9961 lEischer et'all 
l2000h . The measurements fro m the Sloan Digital Sky Sur¬ 
vey (SDSS; lYork et akl 120001) . binned as a function of the 






































2 


ASSEMBLY OF CENTRAL GALAXIES 


galaxy luminosity (iMcKav et al.ll200ll) . allow one to extract 
the masses of the host halos of t he foreground gala xies, simi¬ 
lar to lGuzik & Sell aid ( |2002^ and lYang et al.l ( l2003h . 

In galaxy-galaxy lensing, one measures tangential shear 
around foreground galaxies as a function of the project dis¬ 
tance R. The tangential shear is the ratio of the excess surface 
density to the critical surface density, 7 t(R) = AI](R)/Sciit- 
The critical density equals Ecut = c^T>s/47rGDiDis, where Ds 
is the distance to the source, Di is the distance to the lens, 
and Dis is the distance from the lens to the source. The ex¬ 
cess is dehned as the difference between the local and the 
average surface density, AE(r) = E(r) - E, where E(R) = 
2R~^ Jq T,(R')R'dR' is the mean surface density within R. 

The surface density can be expressed in terms of 
the galaxy-mass correlation function, fgm(r), via E(R) = 
p f ^gmVR^+z^dz- The correlation function between galax¬ 
ies and dark matter can be related to the cross power 
spectrum between the two through a Fourier transform, 
£(r) = (27r)~^ f P^me'^^ d^k. Within the halo paradigm 
dCoorav & Shethl2002^ . the cross power spectru m can be con¬ 
struct ed from the halo mass function dn/dM (iJenkins et alJ 
120011) . the density prohle of a halo p(r,M), and information 
on how the halos are distributed with respect to the linear den¬ 
sity fluctuations. We assume that the density profile is given 
by the NEW functi on dNavarro, Frenk, & Whiteil997t but see 
iMerritt et alJ2005l) . The concentratio n parameter of the N EW 
profile is a function of the halo mass dBullock et alJ200ll) . 

To characterize multiple galaxies sharing the same halo, 
we specify number of galaxies of a given luminosity in a 
halo of a given mass, which is the conditional halo occu¬ 
pation number. We employ the most recent descript ion of 
the halo occupation disttibution dKravtsov et al.ll2()()^ . We 
distinguish between central and satellite galaxies. The satel¬ 
lites galaxies are described by a mean occupation, (Asat) = 
[(M-Mo)/Mi]^, where the parameters Mq ~ ( 10 ^^- 1 O'^)M 0 , 
Mr ^ and A ~ 1, are given in Table 1 of 

IZheng et alJ d2004l) . (5n the scales of interest, the occupation 
number of central galaxies is unity, (Aren) 1- 
The contributions to the cross power spectrum from central 
and satellite galaxies are, 

1 f M dn , , 

Egm.cenW=- / Aeen M(k,M)ufM, 

ric J p dM 

^ ^ ^ {Ns,Mk,M)ug(k,M)dM. (1) 

n^JpdM 

Here, u(k,M) and Ug(k,M) are the normalized Fourier trans¬ 
forms of the halo and galaxy density profiles, respectively. 
For simplicity, we set Ug{k,M) = u(k,M). We also define 
the average density of mass, p = J{dn / dM)MdM, and galax¬ 
ies, ill = J{Ni)(dn/dM)dM, where i distinguishes central and 
satellite galaxies. We have ignored the large scale clustering 
between halos (e.g., there can be a correlation between galax¬ 
ies in one halo and dark matter in another). This term is not 
a significant contributor to the galaxy-mass correlation at the 
radii of in terest (< l/i~* Mpc) as it changes the mass estim ate 
by < 1% dGiizik & Seliakl20 04jMaridelhjnimje^ 

In Figure fTh. we show the lMcKav et aiJ hnml) measure¬ 
ment of excess surface density, separated into four z'-band 
luminosity bins. Overlaid are the best fit models for AE(R) 
calculated as described above. Given that the concentration 
parameter is a function of mass, the fitting is done over a sin¬ 
gle parameter, the halo mass. We concentrate on the z'-band 


data because the relation we expect to derive between the cen¬ 
tral galaxy luminosity and the halo mass is studied at galaxy 
cluster scales in the A-band. 

In Figure [H), we compile our mass estimates, as well as 
the data from ll.in & Mohrl ii2nnl for central galaxy luminos¬ 
ity in gr oups and clusters, and the tot al luminosities in these 
systems (iLin. Mohr. & Stanfordl2004l) . Logarithmic slope of 
the relation between the halo mass and the central galaxy lu¬ 
minosity is different at low and high mass scales. At the 
low mass end, the fits to the gala xy-galaxy lensing data sug- 
gestLf ocM°'^~°'^ , consistent with lGuzik & Seli^ li2002h and 
lYang et ^il2003l) . At the high mass end, the relation is shal¬ 
lower, Lc (X M° ilFin & Mohil2flfl4) . Such a relation is also 
required whe n interpreting the statistics of the 2dF galaxy 
group catalog JYang et alJ2005ft . 

The relation between the halo mass and the total luminos¬ 
ity of galaxies in the halo, howeyer, remains a power-law oyer 
three decades in mass Aot oc with f3 ~ 0.85. Such a power- 
law is consistent with the dependence of the stellar mass-to- 
light ratio on mass, M/L ^ M'/® dBender, Burstein. & Faberi 
Il992h . At the high-mass end, the power-law could also be 
related to the mass function of subhalos within a halo of 
mass M. With the scaling Ns(m\M)dm ^ trfM‘‘~^dm for 
the number of subhalo s with masses between m and m + dm 
JVale & Ostrikeil2004ft . the total luminosity of the parent halo 
scales as Liot{M) « f Ns{m\M)L^{m)dm oc M"“', where L^{m) 
luminosity of the galaxy in the subhalo, which yanishes in 
subhalos too small to contain galaxies. The slope of the Aot- 
M relation is then consistent with the numerical estimates 
aft! 1.8- 1.9 dOe Fiicia et al.l2003) . 

3. CENTRAL GALAXY LUMINOSITY-HALO MASS RELATION 


We now calculate the rate at which the luminosity of the 
central galaxy grows due to the accretion of satellite galax¬ 
ies in the same dark matter halo. The accretion is moderated 
by the dynamical friction time scale, fdf, on which a satellite 
galaxy sinks in the potential of the primary halo. The effec- 
tiye total mass of a satellite galaxy is augmented by the bound 
dark matter mass left behind from the epoch when the satel¬ 
lite was at the center of an isolated halo. The central galaxy 
luminosity thus grows at the rate 


dL 

dt 


dN A 


-t/Ts, 


( 2 ) 


/q dR^ tdf(Ts) 

where N{L^) is the number of satellite galaxies with luminosi¬ 
ties less than L, (the integrated luminosity function). 

A satellite halo experiences torque Tdf = |r x Fdf|, 
where Fdf is the force o f dynamical friction, |Fdf| = 
47rfrifG^Ms(rG^p(r)ln(A)/v? dChandrasekhad 1194^ but see 
iTremaine & Weinberglll984l) . For circular orbits, the yeloc- 
ity of the satellite is the circular yelocity in the primary halo 


Vs = •\/GM(r)/r; M(r) is the mass of the primary halo con¬ 
tained within this radius; p(r) is the density; Ms(rs) is the 
mass of the satellite contained within the radius rj, measured 
from the center of the satellite, at which the satellite is tidally 
truncated, and In A is the Coulom b logarithm. In numerical 
simulations of satelli tes in halos d Velazquez & Whit3ll999t 
iFellhauer et al .1120001) . I n A ft 2; for an explanation of small 
In A see Appendix A of iMilosavlievic Ac Merrittl d2flfl1h . The 
numerical factor /df is of order unity and depends on the or¬ 
bit of the satellite and on the orbital phase space distribution 
of dark matter. The satellite is tidally truncated at the radius 
where its ayerage density equals that of the host halo within 
the orbit of the satellite, Ms{r^)/rl. 
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Fig. 1. — (a) The galaxy-mass con'elation functi on, expressed as the excess surface density AE(/?) as a function of the projected radius R. We show the SDSS 
z'-band data binned into four luminosity bins (from lMcKav et alJ2Q01D . and fits based on the halo model, as described in §[J The model fitting procedure allows 
us to estimate the average halo mass in which the foreground galaxies are contained in each of the the four luminosity bins, (b) Central galaxy luminosity as 
a function of the halo mass. The data points at the low mass en d {large circles) w ere determined by model fits to the SDSS galaxy-mass correlation function 
(see §[3- For comparison, we also show the masses estimated bv lYang et alJ I2Q5^ using the same galaxy-mass c orrelation function {squares), based on direct 
NFW profile fitting to the data or on the scaling of the mass within 260 kpc, as determined bv lMcKav et all 120011) . to the virial radius obtained from numerical 
simulations. At the high mass end, we show a direct measurement of galaxy luminosity and halo mass for a sample of galaxy groups and clusters from Lin & 
Mohr (2004; small circles). We also show the total luminosity of galaxy groups and clusters based on data of Lin et al. (2004; diamonds). The dot-dashed line 
is a power law fit to the combined lensing and Lin et al. data with logarithmic slope (3 = 0.85. The curves, from top to bottom, our models for the central galaxy 
luminosity-halo mass relation with the parameter /i = (0.05.0.1,0.2), respectively, as discussed in §□ 


Let 7s = \/GM{r)r be the specific angular momentum of the 
satellite. The satellite spirals toward the center of the primary 
halo on a time scale fjf = {Ms,/Tif){dJs/dlnr). In the outer 
parts of halos described by the NFW profile, dlnJ^/dlnr « 
0.7- 1. The time scale is the longest when the satellite has 
just entered the virial radius of the primary halo and thus we 
evaluate fdf at r = r^ii, the average densities of the two halos 
when they are touching at the virial radii, p, are both equal 
to 200 times the universal matter density, while the density 
at radius is a factor of ^ (4-7) smaller than the average 
density, depending on the halo concentration. We denote the 
ratio p//o(rvir) = p « 5. If the orbit of the satellite is eccentric 
when it first enters the primary (a very likely circumstance!), 
the effective density that gives rise to the dynamical drag is 
larger than at the edge of the halo. We do not explore these 
complications and subsume the dependence on the halo orbit 
in the factor fa. The dynamical friction then equals 

, _ v _ 

2\/3tt /dfln a ,/GpMfi 


Therefore, fdf exceeds the dynamical time of the primary halo 
^ {GpY^!^ by a factor proportional to the mass ratio of the 
two halos. 

The luminosity function of galaxies sharing the same dark 
matter halo can be modeled by 


dN. 
dL 


(A) = $* 



(4) 


where L* 8.3 x 1 0is a characteristic lumin osity scale, 
-1.3 ^ a < -0.8 (iTin. Mohr. & .Stanfordl 1200^ and refer¬ 
ences therein; a -1 is appropriate for field galaxies and 
a ~ -0.8 for clusters), and is a normalization factor set by 
the requirement that the total luminosity in the cluster equals 


Ltot{M). Thus, $* = Ltot(M)/L^(2-l-a). Substituting equa¬ 
tions m and 0 into equation (|2j we obtain 

dL _ 2\/3^/dfln A Ltot(M) y/Gp 
dt T{2 + a)ri Ll M 

X / e-^^I^M,{L,)L,dL,, (5) 

where M^{L^) is the parent dark matter mass associated with 
the satellite of luminosity L,, and L(M^Y} Y A(M) is the lu¬ 
minosity associated with the smallest satellite, having mass 
Mmin(Af), that can be accreted to the halo ce nter in a growth 
time o f the primary halo, 1 / tHq, where t ^ 1 dWechsler et al.l 
unni- The minimum mass can be then be obtained from the 
condition fdf(Ms = We divide equation Q by 

dM/dt = tMHq to obtain the rate of increase in luminosity of 
the central galaxy per increase in the halo mass as 

dL _ 30/dflnA Ltot(M) ^ 
dM T{2 + Q)riT Ll 

X / (6) 

where we have used the virial density relation to critical den¬ 
sity p = 200 X I^i-kG, where « 0.27. 

To obtain a rough idea about the meaning of equation 
we note that we are interested in the L-M relation for galaxies 
brighter than L*. Then the integral can be approximated by 
the integrand evaluated at the lower limit and multiplied by 
L*. Since = pM with p = pr/SOs/n/Zf^flnA, 


dL 

dM 


1 Ttot(M) 

r(2-l-a) M 


L(pM) 

L* 


l+a 




(7) 
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4. RESULTS AND DISCUSSION 

We integrate the delay differential equation Q numeri¬ 
cally from small to large masses, where at the small masses 
we assume that L{M) = Ltot(M), where LUM) = U{M/M,f 
with M* « 3.5 X IO'^Mq. In the regime L(fiM) <C L* where 
the exponential factor can be neglected, we obtain L{M) ~ 
i.e., the luminosity of the central galaxy is a 
power law, similar to that describing the total luminosity of 
galaxies in the halo. On the scale M^it ~ M*//r the power- 
law behavior breaks, and the growth of the luminosity is sup¬ 
pressed. 

Using a = —l, (3 = 0.85, t = 1, p = 5, Om = 0.3, and In A = 2, 
we estimate /i « 0.16. As discussed above, the evaluation of 
the exact value of ^ is beyond the scope of this analysis and 
may vary from one halo to another. We compare the predic¬ 
tions of the model by plotting L(M) obtained by direct inte¬ 
gration of equation m over the data in Figure^. Because of 
the uncertainty in the precise value of /i, we present curves for 
/r = (0.05,0.1,0.2). For p = 0.1, the critical mass scale of the 
luminosity growth retardation is M^it « 3 x IO'^Mq, which is 
naturally associated with the halos in which the time scale on 
which the satellites merge with the central galaxy is equal to 
the age of the primary halo. This is an additional fundamental 
scale characterizing galaxy formation. This scale is unlikely 
to be directly observed in the luminosity function as the sta¬ 
tistical averaging over several decades in halo mass, implicit 
in the evaluation of the function, erases its signatures. 

As evident in Figure [Qt, Mcut, identihed with the break in 
the L-M relation, is to some degree sensitive to the value of 
/r, which is in turn sensitive to the precise calibration of the 
dynamical friction force, as well as to the orbit of the satellite 
at the point of initial entry into the primary halo. The value of 
/r is also sensitive to r, the dynamical age of the primary halo. 
Both uncertainties plausibly give rise to a factor of 2 variation 
in /i. This can explain the large scatter in the observed cen¬ 
tral galaxy luminosity-halo mass relation in luminous groups 
and clusters. Except for a few outliers, the derived relation 
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reproduces the data over three decades in mass. 

The success in deriving the relation from first principles and 
without reliance on detailed numerical simulations or semi- 
analytical models encourages us to propose several general 
conclusions. In semi-analytic models of galaxy formation that 
are based on efficient gas cooling, the high mass end of the lu¬ 
minosity functi on is generally over abundant (“the overcooling 
problem;” e.g.. lBenson et alJ2003l) . The same high-end of the 
luminosity function is dominated by central galaxies above 
L*. Given that we were able to describe the luminosity growth 
of such galaxies through dissipationless merging, it is unlikely 
that additional cooling inside groups and clusters contributes 
significantly to their stellar mass budget. Thus, gas cool¬ 
ing must be generally suppressed, either through feedback or 
heating, over a wide range of mass scales ranging from clus - 
ters down to L* galaxies (see also iMaller & Bullocldl2004l> . 
While the exact mechanism of cooling suppression remains 
a mystery, we believe that it is not restricte d to clusters w ith 
temperature above ~ 1 keV, as suggested bv iFahi aiil (l2004ft . 

Our results support the paradigm for the formation of gi¬ 
ant galaxies with two fundamental scales, namely, that re- 
lated to the efficiency of shock-heating ^ (1-6) x IO^Mq 
( fDeker& Birnboiriill2004ll . and that related to the efficiency 
of merging ~ (1 -6) x IO'^Mq. The standard M* galaxies, 
which correspond to L* galaxies on the luminosity scale, be¬ 
long between these scales. Since M^it, interpreted as the scale 
at which the L-M relation exhibits a break, is a consequence 
of the drop in the merging efficiency, it also may be possible 
to explain and L* on the basis of dissipationless merging. 
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